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' We present the results of a detailed analysis of a general, unstructured adiabatic quantum search 

, of a data base of A'' items. In particular we examine the effects on the computation time of adding 

' energy to the system. We find that by increasing the lowest eigenvalue of the time dependent 

OA ' Hamiltonian temporarily to a maximum of cx v^, it is possible to do the calculation in constant 

^ , time. This leads us to derive the general theorem which provides the adiabatic analogue of the \/iV 

• bound of conventional quantum searches. The result suggests that the action associated with the 

' oracle term in the time dependent Hamiltonian is a direct measure of the resources required by the 
adiabatic quantum search. 
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I. INTRODUCTION 



in 
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' The subject of quantum computation and quantum information theory has attracted a great deal of attention in 
' the recent past (for reviews, see 0, 0; 0)- This is mainly due to the existence of many algorithms showing that 
' the principles of quantum mechanics can be used to greatly enhance the efficiency of solving specific computational 
I problems. Notable among them are the factorization algorithm due to Shor 0, the algorithms of Deutsch and Jozsa 
1^15 and the data search algorithm due to Grover j^, which is the subject of the present paper. 

Assume that there are a total of N items in a completely unstructured database, out of which exactly one is marked 
Oh, (item m say). The task is to find that marked element in as few steps as possible. Classically, on average N/2 
^ ' steps are required, but Grover was able to devise a quantum search algorithm which finds the marked element in 
^ . 0{^/N) steps [21] . Moreover, this is the best that one can do since Grover's algorithm saturates the lower bound found 
by Bennett et al Q (henceforth referred to as the BBBV bound. See also |9|, llO|, ) . Recently it was suggested 
0, 01 that such gates can be replaced entirely by a continuously time varying Hamiltonian which evolves a chosen 
initial quantum mechanical state directly to the required marked item/state in an efficient way. The details of this 
"adiabatic" approach are as follows: consider TV = 2" items in the database, each associated with one vector in the 
complete orthonormal vector set — 0, ...,N — 1} in the Hilbert space of n spin-1/2 objects. Assume that the 



' initial state of the quantum computer is the symmetric (normalized) state: 



N-l 



N-l 
1=0 

Let I to) denote the eigenstate associated with the marked item. Now define the two Hermitian operators 

Ha = /- |V'o)(^o| (2) 
Hi = I - \m){m\ (3) 

whose ground states are l-^o) and |to), respectively, and the time dependent Hamiltonian 

H{t) = {\- s{i))Ho + s{t)Hi . (4) 
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In the above, s(t) is some monotonically increasing function of time t, subject to the conditions 

s(0) = , s{T) = l . (5) 

where T is the total computation time. Note that Hi, which projects out the marked state, plays the role of the 
"oracle" of the Grover algorithm. 

As long as the time evolution is slow enough, the adiabatic theorem guarantees that the state jV'o) will evolve into 
the desired state \m) in time T. In Ref.|l.ll| it was shown that the optimal running time for such an adiabatic quantum 
search goes as 0{^/N), consistent with the BBBV bound. 

The purpose of this letter is to present the results of a detailed and rigorous analysis of the general form of this 
search algorithm. The corresponding equations reduce exactly to a two-dimensional system which can easily be solved 
numerically for arbitrary N without making any approximations. In principle, the generalized algorithm can be used 
to speed-up the search, but, as one might expect this comes at a cost: it is necessary to inject a large amount of 
energy, albeit temporarily into the system. We illustrate this mechanism with a specific example, and then derive the 
adiabatic analogue of the BBBV bound and comment on its implications. 

II. A GENERAL CLASS OF MODELS 

We begin by considering a time dependent Hamiltonian of the following general form: 

H{t)^fit)Ho+g{t)Hi , (6) 

where Hq and Hi are given by Q and f(t) and git) are arbitrary functions of time, subject to the boundary 
conditions: 



/(O) = 1 
and f{T)^0 



5(0) - 
9{T) = 1 



(7) 

(8) 



Once again, the system evolves from the initial state \ipo) to the marked state |m) in time T. Note, however, that 
unlike the usual algorithm, / and g are not necessarily monotonic (and hence bounded above by one), nor do they 
not obey the constraint f + g = 1 . This will play an important role later on. 

It is a straightforward exercise to find all N eigenvalues {Ej{t)} of the Hamiltonian H{t) and their corresponding 
eigenvectors | Ej , t) . The two lowest eigenvalues are 



(f + g)±^l{f-g)2 + -fg 



The higher energy eigenvalue is (N — 2) fold degenerate: 

m)=f+g 



(9) 



(10) 



From the expression (|5J| it follows that the eigenvalues E± are non-degenerate as long as / and g are non-negative 
and do not vanish simultaneously. 

Without loss of generality we now write the solutions \ip(t)) to the time dependent Schrodinger equation in terms 
of the complete basis of energy eigenstates: 



\m)= E '^jit)e'^-fo^'^'>''\E,,t) 



(11) 



The probability amplitudes a±{t) for the two lowest energy eigenstates decouple from those of the higher excitations. 
They obey the simple equations: 



In the above ' = d/dt and 



F. ^ -F*_ ^ 



■ exp I 



jdt' 



(12) 



(13) 



3 



and we have defined: 



io{t) := E+it) - E_ it) ^\lif- gf + -fg 



We have also used the following matrix elements: 



E+,t 



E±,t 



dH 






E^,t 


dt 


dH 




dt 


Ei,t 



VN~T fg-gf 
N uj 

, . 



According to the adiabatic theorem as used in ,12!|, if there is a small number e << 1, such that for all t, 



E+,t 



dH 



dt 



E^,t)— <e 



(14) 

(15) 
(16) 

(17) 



0, then the probability of transition to the excited 



state, |a4_(t)p, will not exceed a number of order e^. Note that for the theorem to hold, the eigenvalues E±{t) must 



and if the system starts in the ground state |-0o) of H{t) at t 

never cross. For a given desired computational accuracy e, Ea. (|17|l limits the rate at which the Hamiltonian can 
evolve, and therefore puts a lower bound on the total running time T. In order to calculate the bound, one needs 
more information about the form of functions f{t) and g{t). For example, if f{t) = 1 — s{t) and g{t) — s{t), as in 
Eq.Q, the total running time T can be minimized by choosing ds/dt to saturate the bound 1)17(1 at each instant in 
time. This results in a running time of order \/N consistent with the BBBV bound Q. 



III. A SPECIFIC EXAMPLE 



We now illustrate with a specific example that more general choices for the functions f{t) and g{t) allow the 
computation to be performed considerably more rapidly. Consider f{t) and g{t) to be quadratic functions of s(t): 

f{s) = I - s + As{l^ s) , g{s) ^ s + As{l - s) (0 < s < 1) , (18) 

where A is a constant and s{t) is again a function of time that varies monotonically between and 1, so that the 
required boundary conditions (|7I8|I are satisfied. Note that / and g are not monotonic in this case. Their peak values 
are A/4+ 1/2 =F l/AA respectively (occurring at s = 1/2(1 =F l/A) respectively). The "adiabaticity condition" p7|) in 
this case becomes (we use d/dt — (ds/dt) d/ds) 



dsy/N^ l + A(l-2s + 2s2) 

<^ ^ (19) 

^ [(1 - 2s) + ^s(l - s) [1 + A - s(l - s)A^]\ ^'"^ ~ 

For a system satisfying the constraint 1)19(1 . the minimum computation time can again be found by adjusting s(t) 
so that the inequality is saturated at all times. This yields a minimum running time of: 

I VF^ |l + A(l-2,s + 2s^)| 

Tmin = / ds -TT- . (20) 

« ^ h [(l-2s) + ^.5(l-s)[l + A-s(l-s)A2]]^/^ 

If A = 0, the result of 12] is reproduced. However, by choosing A appropriately, considerably lower running times 
can be achieved. For example, A = ypN 3> 1 we get : 



N 

The limit 1(21(1 can be obtained by Taylor expanding the integrand in ((20(1 about N = oo. This is possible for the 
choice A = ^/N because the denominator does not vanish anywhere when N = oo. One has to be considerably more 
careful in deriving the limit for other values of A. We have also verified 1(21(1 numerically. 

Thus we see that the data search can be completed in a constant time, independent of N. The cost of this 
improvement is that the energy of the system gets very large for a finite time before returning to zero at the end of 
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the computation. In particular, for the choices of / and g in H18|l . E-{s) reaches a maximum of {A-\- 1)/4 (for N 1) 
at 5 = 1/2, so that \i A — -s/ZV, the maximum energy grows with ^/N as well. 

By numerically solving (|12|1 . it can be confirmed that the quantum computer indeed remains in the ground state 
at all times. One can explicitly calculate the transition probability P-{t) = |a_(t)p as a function of time (by the 
equations of motion (|12|1 . the corresponding -P+(t) = |a^(t)p is guaranteed to be 1 — P-{t)). Figure illustrates a 
sample solution with A = ^/N and s{t) chosen to saturate the inequality (|19|) at all times t. The transition probability 
|a_p remains above ^ (1 — e^) and the graph is plotted till the end of running time = 1 + 7r/2 ^ 2.57. 

Time dependence of P_(t) 
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FIG. 1: Numerical solution for P_(t) = |a_(t)p for = 100, 000 and e = 0.002 with A = y/N. 



We have explored numerically other choices of A. For A ^ A^", where a > 0.5, the calculational running time 
approaches a constant for large N, with the constant of proportionality decreasing as a increases. A systematic 
analysis of this problem is in progress , but in light of this and the discussion below, it is interesting to speculate 
what happens to the running time for a non-polynomial choice such as ^4 ~ e^. 



IV. ADIABATIC VERSION OF THE BBBV BOUND 



At this stage it is natural to ask what are the most general conditions that allow a speed-up of the search algorithm. 
This question is answered in the following theorem, which generahzes the results of [M 1^. ITol ITU IT^ . 

THEOREM: Suppose there exists an unstructured adiabatic quantum search algorithm that can success- 
fully evolve the initial state ^ to any one of a complete set of orthogonal states via a Hamiltonian 
containing an oracle term of the form g{t)\m >< m\. Then the time T for the computation is bounded 
below by the relationship: 

1 /"^ . > , kv^ 



^ , 9it)dt > (22) 

in the limit that N is large {k is a constant of order unity). 

This inequality is virtually identical to the one obtained by lll^ . The one important difference is that we do not 
restrict g{t) to be less than one, nor do we impose the constraint f{t) + g{t) = 1. Details of the proof can be found 
in |l3j| so we will not repeat the details, merely outline the important ingredients. The key is to decompose the 
Hamiltonian into two parts: 

H{t) ^ H2mit) + Hi{t) , (23) 

where 

H2m{t) = -g{t)\m){m\ . (24) 
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contains all the dependence on the needle state |m) and Hi{t) is an arbitrary term independent of the state |m) (in 
particular, for the Hamiltonian ©, Hi{t) = f{t) + g{t) — /(t)|V'o)('0o|)- Note, however, that Hi{t) is completely 
general, apart from the restriction that it not depend explicitly on the state |m). 

Consider two computers at time t, evolving to states \m) and |m') respectively, represented by the wavefunctions 
\ipm,t) and \ipm',t) respectively subject to the boundary conditions: 

|^™,0) = |^™',0) = l^o) (25) 
|V'™,T) = |m) ; \i,^,,T) = \m') . (26) 

In terms of the specific model we considered earlier, this requires f{T) — and g{T) = 1. Following the steps in 
it is straightforward to obtain the inequality: 

A]\f3/2 i-T 

^ [i-|(^™,r|V'™',r)n <^— y g{t)dt ■ (27) 

The assumed orthogonality of the final states at time T implies that 

1- |(V'„,T|V'™',T)|2 > fc ,Vm^TO' (28) 

where fc is a number of order 1. Summing over m and m' in Eg. (1271) fo r ^ 1 yields Ea. (|22|l as desired. Again 
we stress the Eq.ljJD) is virtually identical to the result obtained in [l^- The key difference concerns our relaxed 
assumptions about the functions / and g in the Hamiltonian. 

Eq.|23) shows that the BBBV bound can only be beaten if the mean value of g{t) over the time T grows as some 
power of N. If this power is 1/2 then the above theorem states that the running time T may be bounded by a constant 
independent of TV, as in the example quoted in this paper. Note that this result is a generalization of a similar bound 
obtained in for the special case given by Q, in which g{t) never exceeds order unity. It also is a natural extension 
of the result presented by Farhi and Gutmann in 9] for constant g. 



V. CONCLUSIONS 



In summary, we have presented the results of an analysis of a generalized adiabatic quantum search algorithm. 
The corresponding Schrodinger equation was shown to reduce exactly to a two dimensional system for arbitrary 
N. We derived the adiabatic analogue of the BBBV bound. Our theorem shows that the optimal speed normally 
associated with Grover's search algorithm can be improved in this framework by a suitable choice of the time-dependent 
Hamiltonian. As one might expect from dimensional grounds this speed-up requires an increase in energy, at least 
temporarily. However, it should be emphasized that it is not the total Hamiltonian that needs to be increased, only 
the coefficient in front of the oracle term. In principle this leaves open the possibility of speeding up the search while 
keeping the ground state energy of the system small. Another way to keep the ground state energy zero would be to 
use a new Hamiltonian obtained from by subtracting the term E-{t)I from the latter. Note that although the 
resultant Hamiltonian cannot be written in the form JHq + gHi, since £'-(0) = = E-{T), it would still evolve the 
initial state to |m). Moreover, the 'gap' 0(t), as well as the matrix element (|15|) will remain intact, implying that 
the running time will still be a constant, given by (|21|) . In any case, our analysis suggests that the physical quantity 
J- g(t)dt provides an adiabatic analogue of resources required for the unstructured search, just as the number of 
operations does for the conventional quantum search. 

It must be emphasized that the system considered here is highly idealized. The Hamiltonian is non-local in the 
sense that Hq and Hi require all qubits to be coupled simultaneously. It is therefore not clear how to implement 
such a Hamiltonian in a realistic physical system. One should therefore investigate the circumstances under which 
one can find simpler, more local, Hamiltonians that have the same ground states and hence can be used as the basis 
of a realistic adiabatic quantum search algorithm. This issue will be the subject of a separate publication jl6||. 
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